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Abstract 

In steady-state non-isentropic flows of perfect fluids there is always 
thermodynamic generation of vorticity when the difference between the 
product of the temperature with the gradient of the entropy and the gradi- 
ent of total enthalpy is different from zero. We note that this property does 
not hold in general for complex fluids for which the prominent influence 
of the material substructure on the gross motion may cancel the thermo- 
dynamic vorticity. We indicate the explicit condition for this cancellation 
(topological transition from vortex sheet to shear flow) for general com- 
plex fluids described by coarse-grained order parameters and extended 
forms of Ginzburg-Landau energies. As a prominent sample case we treat 
first Korteweg's fluid, used commonly as a model of capillary motion or 
phase transitions characterized by diffused interfaces. Then we discuss 
general complex fluids. We show also that, when the entropy and the to- 
tal enthalpy are constant throughout the flow, vorticity may be generated 
by the inhomogeneous character of the distribution of material substruc- 
tures, and indicate the explicit condition for such a generation. We discuss 
also some aspects of unsteady motion and show that in two-dimensional 
flows of incompressible perfect complex fluids the vorticity is in general 
not conserved, due to a mechanism of transfer of energy between different 
levels. 

1 Introduction and statement of the problem. 

For steady-state non-isentropic flows of perfect fluids, the following pointwise re- 
lation holds throughout the regular region B of the three-dimensional Euclidean 
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point space £ occupied by the fluid: 



u> x v = dgradn — gradH, 



(1) 



where v is the velocity of each material element (v = v(x)), u> the vorticity 
(u> =curlv), $ the temperature, 77 the entropy density, H the total enthalpy 
density (grad and curl are calculated with respect to x GB). Equation (1) is 
known as Crocco-Vaszony's theorem (see [1-3] ) . It states basically that we may 
observe thermodynamical generation of vorticity when 



For complex fluids, for which the gross behavior is influenced by a prominent 
presence of material substructure, the condition (2) is no more sufficient to 
assure vorticity and we may have situations in which the effects due to the 
substructure cancel the thermodynamic generation of vorticity, or may produce 
it even when the entropy and the total enthalpy are constant. 

In what follows we give explicit conditions for the occurrence of such a topo- 
logical transition, due to the complex structure of the fluid, from a vortex sheet 
to a shear flow and vice versa 1 . To fix basic ideas, we consider first the simplest 
case of complex fluid, the one of Korteweg, in which effects associated with the 
gradient of the specific volume are accounted for. Then, we treat a general case 
(following the unified theory in [4]) in which each material element at x, at 
the time t, has an articulated morphology and can be considered as a 'system', 
described by a coarse grained order parameter v, taken here as an element of an 
'abstract' differentiable manifold M. to cover a class of complex fluids as large as 
possible. The possible non-uniform distribution of the order parameters influ- 
ence the formation of vorticity during the flow. Finally, we discuss the unsteady 
case and show that in a two-dimensional flow of incompressible perfect complex 
fluids the vorticity is in general not conserved (as in simple perfect fluids) un- 
less an explicit condition is satisfied. The influence of the material substructure, 
rendering not conserved the vorticity, is not of viscous nature: there is just a 
transfer of energy from the coarse grained (macroscopic) level to the substruc- 
tural one, altering the macroscopic energetic landscape of the vorticity (that 
may generate possible drag reduction in turbulent flows). This phenomenon 
implies the possibility that even in a two-dimensional flow of complex fluids a 
vortex may be destroyed and, possibly, reconnected in two or more subvortices. 

For the sake of simplicity we do not consider applied external fields. 

2 The case of Korteweg's fluid. 

Motivated by the work of van der Walls and previous ideas of Lord Rayleigh and 
Gibbs, to model at coarse scale the effects of interfacial tension in capillary phe- 
nomena, in 1901, Korteweg [5] introduced a model of fluid in which the energy 

1 Such a kind of transition may occur in disparate circumstances, as, e.g., geophysical flows 



dgrad-q — gradH ^ 0. 



(2) 



[18]. 
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depends on the gradient of the density. Experimental works have supported 
this approach [6-8] and the model has been source of analytical investigations 
[9, 10]. A detailed critical review about this topic can be found in [11]. The 
work of Korteweg has been the germ of a general modeling of 'materials of grade 
N' in which higher order deformation gradients up to grade N are involved in 
constitutive equations (see [12] for the general thermodynamic treatment). In 
1985 Capriz [13] proved clearly that the constitutive dependence on higher order 
'deformation' gradients is an indicator of the prominent presence of a complex 
material substructure which needs the introduction of some order parameter, 
under the assumption of the existence of an internal constraint linking the order 
parameter with the deformation (Korteweg's model follows as a special case). 
Our prominent interest in the present setting is not properly on capillary phe- 
nomena, the original motivation of Korteweg's model, rather it is focused on the 
circumstance that such a model is currently used to represent phase transitions 
on the basis of a diffuse interface approach (see [14] for a detailed review). In its 
simplest form, in fact, the free energy is of Cahn-Hilliard type: the gradient of 
the mass density account for the free energy excess of the interfacial region be- 
tween phases. Such an excess of energy may influence the generation of vorticity 
in non-isentropic flows: it may cancel or produce vorticity. 

2.1 The main result about Korteweg's fluids. 

Here we analyze the simplest case of complex fluid, i.e., the perfect Korteweg's 
one: a compressible perfect fluid characterized by the functional dependence of 
the potential <f> on the gradient grade of the specific volume t (being i = I (x) = 
detg (x), with g (x) the metric in space at x 6 B, and grad is calculated with 
respect to x) in addition to I itself; in other words, we consider the case in which 
(f> = <f> (i, gradi, rf), with (f> a sufficiently smooth function. 

In the following, we write formally <f> (i, gradi, rf) being conscious that since 
4> is a scalar function depending on a vector (namely gradi), for reasons of 
objectivity (frame indifference under the action of 5*0(3)) <j> must depend on 
the vector gradi only through its amplitude \gradi\. The simplest choice of 4> is 
Cahn-Hilliard's form: <j> = f (i, rf) + hf3 \gradi\ , with /3 some constant and / a 
function describing, e.g., a two- well energy; then <f> may account for the spinodal 
decomposition. We also indicate with p the current density of mass of the fluid 
and consider the referential density p equal to 1 so that i = p^ 1 . In general, 
the Lagrangian density of the Korteweg's fluid is considered to be given by 

1 2 

7jPl v l + X - P4> \h gradi, rf) (3) 

for non-isentropic flows, where x is a sufficiently smooth function such that 
®iiX 7^ 0, x ( L i — and the equality sign holds when i — 0. The function 
X accounts for possible substructural inertia effects that may generate gradient 
inertia. In other words, the appearance of the time derivative of i, namely i, 
in @ is not of viscous nature, rather it is of inertial nature and the physical 
dimensions of x are the ones of a kinetic energy. Notice that the kinetic energy 
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is the partial Lcgcndre transform with respect to i of x which is called kinetic 
co-energy, following the terminology of [4]. The possible existence of gradient 
inertia effects in the Korteweg's fluid has been pointed out in [13] , however, these 
effects seem to be negligible unless the fluid oscillates at very high frequencies. 
In any case we formally account for the presence of x m the following, because 
in the general setting treated below its effects may be sensible. 

Hereafter, the dot indicates the scalar product associated with the natural 
dual pairing; moreover, we denote with d y the partial derivative with respect to 
the argument y and, for different types of linear operators, with a superscript T 
we indicate the transposition rendered explicit by a subsequent detailed expres- 
sion in components, when the meaning of the transposition is not immediately 
evident. For A and B tensors of the same order, A ■ B indicates scalar product. 

The main result of the present section is summarized in the proposition and 
the corollaries below. 

Proposition 1 For a steady-state non-isentropic flow of a Korteweg's fluid, the 
following pointwise relation holds: 

oj x v = "Dgradv — gradh — grad (L 2 divpd gra dL4> + d gr adu4> ' gradi) + 

+igrad (JL (d iX ) - . (4) 

In equation h is the total enthalpy density of the Korteweg's fluid and 
its explicit expression is given in Section 2.3. 

Corollary 1. In a perfect Korteweg's fluid, the gradient effects cancel the 
thermodynamic generation of vorticity when 

■dgradn — gradh = grad {i 2 divpd gra di<j> + d gra di,4> ' gradi) — 



-igrad (dix) - d.x) ■ (5) 

Corollary 2. In the case in which both the specific entropy and the to- 
tal enthalpy are constant throughout the flow of a perfect Korteweg 's fluid, the 
gradient effects may generate vorticity: 

u> x v = -grad {i 2 divpd grad i4> + d gra d L <t) • grade) + igrad (d L x) ~ d L x^j ■ 

(6) 



2.2 Preliminaries to perfect Korteweg's fluids. 
2.2.1 Equilibrium. 

In absence of body forces, the equilibrium of the fluid in B can be determined 
by evaluating the variation of the functional 

$ = / p<t)(i,gradL) (7) 
Jb 
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and putting it equal to zero, i.e. 

4> (i, gradi) = 0. (8) 



for any choice of <5x, where 8 indicates the total variation (in the sense used in 
elementary treatises when one deals with 'total' time derivatives), so that we 
have (see [4], p. 31) 



i pcf) = 18 {pip) + / p0n-(5x 
Jb jb Job 



dB 

(8 (p</>) + div Go#x)) = / pU, (9) 
Jb 

where n is the outward unit normal to the boundary dB of B. 
We note that 

8(f> (t, gradi) — d L (f>8i + d^adi^S grade, (10) 

where Sl and Sgradi are the total variations of the relevant arguments calculated 
taking into account that l varies in space. We have in fact 

8b = idivSx., (11) 

as a consequence of Euler formula, and by chain rule 

Sgradi — gradSi — (gradSx) T gradi (12) 

(see [4] for other details). By developing the variation in ©, taking into ac- 
count (|TU|l - ltTT|l and using repeatedly the Gauss theorem, we recognize that the 
arbitrariness of <5x implies in the bulk the pointwise balance 

divT = 0, (13) 

(as Euler-Lagrange equation) with 

T = -pi - T E , (14) 

where I is the second-order unit tensor, p a non-standard pressure given by 

p = -pid L 4> + idivpdgr a di4>, (15) 

and T E a non-viscous stress of Ericksen's type called in this context Korteweg 
stress, i.e. 

T E = gradi <g) pd gracU (l>- (16) 
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2.2.2 Flows. 

Let us consider first isentropic flows. The dynamics of a Korteweg's fluid can 
be simply derived by making use of a Hamilton variational principle of the form 

sj J Q|v| 2 +xM)) -J 5 J/xl>(L,grcuk)=0, (17) 

where the time t runs in [0, t] and 5 means total variation as before. By devel- 
oping i|17|) . we realize that the arbitrariness of <5x implies the standard balance 

p^- = divT. (18) 

The stress tensor T is of the form 

T = —pi — T E , (19) 
but now the non-standard pressure p accounts for kinetic pressure effects and is 

p = -pid L 4> + idivpd graded - P L ~kdiX + d tX, (20) 
while T E remains the same of <|16[) . 

2.3 Proof of Proposition 1. 

We denote with q the amplitude of the velocity, i.e., q = |v|. Then, by using 
the standard Lagrange formula 

^ = <9 t v+ {gradw) v = <9 t v + u> x v+^gradq 2 , (21) 

we recognize that in the case of steady motion, i.e., when dtv = 0, the balance 
equation (|18f) can be written as 

u> x v = ——gradq 2 — igradp — idivT E . (22) 

When the flow is not isentropic, we may identify the potential 4> with the 
internal energy and write 

<p = <p(i,gradL,r)), (23) 

being i] the entropy density. From l|23() . we define the specific enthalpy density 
£ as the opposite of the partial Legendre transform of <p with respect to the 
kinematic variables, so that we have 

(, = ()>- id L (f> - d gradL (f> ■ grade, (24) 

i.e., from lf2U|l. 

£ = (ft + ip - i 2 divpdgradi4> + Lp- dg r adi4> • gradi, (25) 
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where 

P = Pt— - d t ,x (26) 

is the kinetic pressure due to possible gradient inertia effects. By indicating 
with ■& the temperature given by d n (j), the calculation of the gradient of Q25|l 
allows us to obtain 

— igradp = 'dgrad-q — grad£ + igradp — L 2 grad {divpd gra dL4>) — 

- {graddgradi,4>) T grade - t (divpd grad L<j)) grade. (27) 

By inserting l|27|) in 1|22|) . taking into account that 

divT E = (divpdg r adt(l>) gradi + (pd gra dL4>) grad 2 L, (28) 

and defining the total specific enthalpy h of the Korteweg's fluid as 

h = \q 2 +^ (29) 

we obtain and Proposition 1 is proven. The corollaries are immediate con- 
sequences. 

3 General complex fluids. 

We follow here the unified setting proposed in [4] (and further discussed in [15] 
adding the treatment of singular surfaces) with the aim to encompass a class 
of complex flows as large as possible. We then consider each material element 
as a subsystem collapsed into a place x and described by a coarse grained 
order parameter v. We do not specify v; we require only that it belongs to a 
differentiable paracompact manifold M. without boundary where we presume 
that physical circumstances induce a single choice of metric and of connection. 

We have then a sufficiently smooth mapping B 3 xh^v = i> (x) G.M. 

We account for gradient effects, i.e., we consider the interaction between 
neighboring material elements under the suggestions of the results about Ko- 
rteweg's fluid. These weakly non-local effects are the crucial mechanism altering 
the standard thermomechanical behavior of the vorticity. 

Under previous assumptions, the Lagrangian density of a general complex 
fluid can be written as 

1 2 

-p|v| + x (y, i>) - P4> (i, v, gradu, rj) (30) 

where the substructural kinetic co-energy x (if perceptible in experiments) de- 
pends now on v and its rate, and the potential <fi accounts for weak non-local 
gradient effects. The contemporary presence of u and gradu in the constitutive 
list of entries of (j> is not simply matter of choice of generality in the modeling 
(due, say, to a principle of equipresence) , rather it has geometric origin. The 
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manifold A4 is in general not trivial (in most cases of physical interest it is not 
a linear space); the pair (v, gradv) collects the peculiar elements of the tangent 
mapping Tv : TB — > TM between the tangent bundle of B and the one of 
M.. In general, the two elements cannot be separated invariantly unless A4 is 
endowed with a parallelism. One asks also a physically significant parallelism. 
Even when Ai is Riemannian, there is no immediate physical reason, in prin- 
ciple, prescribing that the consequent Levi-Civita connection over Ai have a 
prevalent role even in exotic cases. Taking into account the possible articu- 
lated geometry of A4, in general we account for both v and gradv as entries of 
<f) when we decide to model interactions between neighboring material elements. 

We list just below some well known special cases of our treatment to clarify 
the topic; other cases are collected at the end of this section. 

1. For nematic liquid crystals, v is identified with a unit vector, say n, in- 
dicating the direction of the prevailing nematic order of stick molecules 
at x; Ai is the unit sphere with the identification of the antipodes and <j> 
Oseen- Frank potential (see e.g. [16]). 

2. For semi-dilute polymeric fluids, the material element is like a box con- 
taining a population of polymer chains described by end-to-end vectors 
r and v is identified at each x with a second order tensor, TZ, precisely 
the second order moment of the distribution of the r's averaged over the 
population considered (see e.g. [17]). 

3. The simplest reasonable form of <j) is of Ginzburg-Landau's type, namely 
<f> = j (l, v, rf) + ^a\\gradv\\ 2 , with a some appropriate constant and 
||<7radz/'|| the norm of gradv (x) in Horn (T X B,T V A4), the space of lin- 
ear maps between the tangent space of B at x and the tangent space of 
M at v = v (x). 

Notice that if we impose the internal constraint Z (x) = v (x) we reduce the 
treatment to the Korteweg's fluid. 

3.1 The main result about complex fluids. 

By taking into account the structure (|3U[) of the Lagrangian density, we may 
formulate the natural generalization of Proposition 1 . 

Proposition 2 For a steady-state non-isentropic flow of a general perfect com- 
plex fluid, the following pointwise relation holds: 



u> x v 




(f>) gradgradv, 



(31) 



where h c is the total enthalpy of the complex fluid defined below. 
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By indicating with a components on some chart over M. and with ... the 
standard spatial components, (13 1|> may be read as 

(w x v) ■ = (gradr]) i - (gradh c ) i - 

grad (^divd gradv <j> - ^ (dox) ~ dvX 

- {graddg r adv$)\ a (gradv)™ + l (gradv)" (divd gr adu4>y aj 
+i {dgradv<t>) 3 a (gradgradv)^ . (32) 

Corollary 3. In a perfect complex fluid, the presence of material substruc- 
ture cancels the thermodynamic generation of vorticity when 

dgradrj — gradh c = (gradd gra dv4>) T gradv — i (gradv) T divdg ra dv<t> + 

+ (^grad [divd gradv 4> - ^ (dox) - duX^j v ~ 

-i{d gra dv<i)) r gradgradv. (33) 

Corollary 4. In the case in which both the specific entropy and the total 
enthalpy are constant throughout the flow of a perfect complex fluid, the presence 
of material substructure may generate vorticity: 

u> x v = — (gradd gra( iv4>) T gradv + l (gradv) T divd gra d v <t> + 
+i {d g radv(t>) T gradgradv - 

grad (^divd gra d v 4> - ^ (dox) - dvX^j J v. (34) 

3.2 Proof of Proposition 2. 

3.2.1 Flows of isentropic perfect complex fluids. 

To obtain the field equations for general complex fluids, we consider first the 
total Lagrangian 



L = j B (\p l v | 2 + X ("> £/ ) ~ ( fc > ^ 9radv) 



(35) 



and set equal to zero the first total variation (in the sense of pUfl ') of its integral 
over the time interval [0, £], namely 

6 [ L = 0. (36) 
Jo 

We consider also that 

Sv =5v + (gradv) 5x, (37) 
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in components 
and that 



(38) 



J pS(f> — J p {d L (j}8i + d v (j) ■ Sv + d gra dv<i> ■ Sgradv^j . (39) 

By using repeatedly Gauss theorem, variations vanishing at the boundary 
of B and the guidelines of Section 2.2.1 (see also [4], p. 30-34), we obtain the 
following balances: 

P^jr = divT , divS z = -j.^uX - 9uX, ( 40 ) 
where the stress tensor T has in this case the form 



T =pid l 4>l - (gradu) d gradu (f> (41) 



and we use S and z as 



S = pdgradui 1 , z =pd„4>, (42) 

calling them microstress and self-interaction, respectively. The pointwise bal- 
ance of substructural interactions (140b ) states that self-interactions z accrue 
within each material element (due to its substructure) to balance the 'contact' 
interactions S between neighboring material elements (due to the relative change 
of substructure). At each x, z is an element of the cotangent space T*A4 while 
S belongs to Horn (T*B,T*A4). As for the pair (i/,gradv), even in this case 
the non-trivial structure of A4 does not allow us to separate invariantly the 
microstress from the self-interaction, in general. However, when an appropri- 
ate parallelism permits us to select an explicit form of cj> depending on the sole 
gradu and substructural inertia effects are negligible, equation (|40b ) reduces to 
a divergence form. For example, such a situation happens in the special case of 
ncmatic liquid crystals when we select just the simplest Frank's potential and 
reduce ourselves to investigate only the physics described by harmonic maps 
taking values in S 2 [19]. 

As anticipated in discussing the special case of Korteweg's fluid, the kinetic 
co-energy x (v, v) is such that its Legendre transform with respect to the rate i> 
coincides with the substructural kinetic energy k (y, 0), i.e., k = doX ^ ~X- It s 
contribution may be in general sensible. Notice that, even when k (f, v) is of 
the form i;> Q fi a ^z>^, the explicit form of \ differs from k by an addendum of the 
type \ a v a , which may be helpful to describe some prominent physical effects. 
For example, in the case of magnetizable fluids, if we select the magnetization 
vector as order parameter and constraint its rate to be of gyroscopic nature, we 
reduce (|4l1b ) to the equation of Gilbert. 
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3.2.2 Details about Proposition 2. 

When the flow is not isentropic, the potential is of the form 

4> = <f){L,v,gradv,ri) (43) 

and we define (as in Q24)0 the specific enthalpy density £ c of a complex fluid as 
the opposite of the partial Legendre transform with respect to the morphological 
descriptors, namely 

£ c = (f> - id L (j) - d v cj> ■ v - d gradv (j) ■ gradv. (44) 

We may also rewrite the term —Ld L <j) as ip, where p is the pressure relevant for 
this case, given by —d u 4>- Equation 1221) becomes 

wxv = — -gradq 2 — igradp — idivT E , (45) 

with 

f E = {graduf d gra d»<t> , (T s )' = (gradu)^ (d grad ^) 3 a . (46) 
Then we get 

divT E — (gradv) T divd g rci dv 4> + {d g radv4>) 7 gradgradv . (47) 
The gradient of (|44|l allows us to write also 

—igradp — dgradrj — grad£ sC — (gradd v (f) T v — (gradd gra dv(j>) T gradv. (48) 

By inserting (|4*8f and iHTf into (f4*5*)l and taking into account fcKh ) and (@2p), 

we get Proposition 2 where 

hc = \q 2 +L- (49) 
The corollaries are immediate. 

3.3 Further special cases 

Well-known cases covered by our general modeling have been listed at the be- 
ginning of Section 3, namely, nematic liquid crystals and polymeric fluids, and 
we may also mention ferrofluids once, in presence of external electric fields, we 
add to <j) the electromagnetic energy and develop appropriate variations of the 
relevant fields. 

We may also investigate more exotic cases by making use of the flexibility 
of the general unifying framework in which we act. Some examples are listed 
below, but they do not exhaust the range of possibilities. 

1. (Smectic-A liquid crystals.) In the smectic-A phase, liquid crystals appear 
organized in layers in which stick molecules tend to be aligned orthogo- 
nally to the layer interface. Within each layer the molecules flow freely; in 
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the orthogonal direction (where the behavior is basically the one of a one- 
dimensional crystal) they may permeate from a layer into another. Linear 
[20] and non-linear [21-23] continuum models have been formulated to de- 
scribe various aspects of the elasticity and hydrodynamics of smectic-A 
phase, paying a few attention on the behavior of possible vortices there. 
Since the current formulation of the hydrodynamics of smectic-A liquid 
crystals falls within the general setting discussed here, under the condi- 
tions of Proposition 2 and subsequent corollaries, we could foresee forma- 
tion or cancellation of vortices. In fact, natural ingredients for modeling 
the smectic-A phase are a vector n £ S 2 at each xe6 describing the 
local orientational order, and a scalar function w parametrizing the lay- 
ers. More precisely, we may write [23] w (x, t, aX), with A an appropriate 
length scale and a running in a set of integers, to account for possible 
unequal spacing of layers. However, at a (gross) coarse grained scale, we 
may imagine that w (x, t, •) be defined in a continuum approximation on 
an interval of the real line so that, as shown in [23], \gradw\~ measures 
the current thickness of the layers. Then, we can select the order param- 
eter as the pair (n, to); however, in absence of tilt (a condition reasonable 
outside the defect core) we have strictly 



gradw 
\gradw\ ' 



(50) 



and the simplest choice for the potential <j> in the compressible case is of 
the form (see [22]) 

1 2 1 2 

4>(l,w, gradw) = <£(i) + (\gradw\ - 1) + --y 2 (divn) , (51) 

with the identification (|50|l . and with j 1 and 7 2 material constants. In 
Q51[). the term (\gradw\ — 1) accounts for the compression of layers while 
(divn) 2 describes the nematic phase and is the first addendum of (three 
constant) Frank's potential (see [3], p. 55). Inertial effects can be con- 
sidered by putting % = 5P a |n| 2 , with a a material constant [23]. The 
special counterparts of (|31|l . I|33|l and l]34p. associated with l|51|l . follow af- 
ter straightforward calculations. To simplify further the matter, one may 
imagine to have an incompressible flow, neglecting 4>(l), a circumstance 
not incompatible with the assumed compressibility of the layers. Taking 
into account the incompressible limit of l|51(l . namely <fi (u>, gradw), the ap- 
plication of Proposition 2, in absence of prominent substructural inertia 
of stick molecules in smectic-A phase, leads to 

u; x v = "dgradrj — gradh c — (gradS) T gradw — 

—wgrad (divS) + udivS (gradw) + lS (gradgradw) , (52) 

with 

S = 7i (\gradw\ — 1) n — 7 2 \gradw\ 1 (I — n ® n) grad (divn) , (53) 
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and I the second-order unit tensor. 

2. (Polyelectrolyte polymers and polymer stars.) In the case of polyelec- 
trolyte polymers suffering possible polarization, we may imagine to assign 
at each point x £ B not only the dyadic tensor 1Z, mentioned at the 
beginning of Section 3, but also a polarization vector p. By indicating 
with Vec the translation space over £, and with B p a ball of Vec cen- 
tered at zero and with radius equal to p, the maximum amplitude of the 
polarization, we have M = Horn (Vec, Vec) x B p and we may construct 
the relevant hydrodynamic theory, once an explicit expression for 4> has 
been selected (see examples in [17]), adding, in presence of external elec- 
tric fields, the electromagnetic energy. When we consider also polymer 
stars, the picture becomes more articulated and we may imagine to have 
A4 = Horn (Vec, Vec) x B p x [0, b]. In this case, at each x £ B, an arbi- 
trary element of the interval [0, b] of the real line describes the radius of 
gyration of the polymeric molecules. 

3. (Superfluid helium 4 He.) The work of Rasetti and Regge [24] allows us 
to manage suitable Poisson brackets to describe 'directly' the dynamics of 
vortices in superfluid helium (see also the developments in [25]). However, 
following the suggestions of [4] (p. 7), we may choose to adopt the general 
point of view used here and to select A4 = S 1 C C, using in this way a 
potential <f> with a Ginzburg-Landau structure with complex entries. The 
application of Proposition 2 and subsequent corollaries follows. 

4 Concluding remarks. 

In the case of Korteweg's fluid, diffuse interfaces modeled by the specific vol- 
ume may generate an energetic wall against the thermodynamic formation of 
vorticity. However, there are conditions for which such an excess of energy may 
favour vorticity even for an isentropic flow. 

Analogous phenomena occur in the case of general complex fluids. Notice 
that, excluding possible inertial terms associated with the gradient of v and 
the derivative of </> with respect to gradv are prominent ingredients into the 
vorticity equation l|31|l . This circumstance suggests that inhomogeneous spatial 
distributions of substructures (preferably with high gradients) in steady-state 
flows may generate energetic barriers to vorticity or may favour it. Basically, 
the energy involved is the one associated with weakly non-local interactions 
between neighboring material elements. 

For non-stationary incompressible two-dimensional flows the vorticity is in 
general not conserved even in the perfect case due to a mechanism of transfer of 
energy between the coarse grained level and the substructural one, ruled by the 
tensor T E . In fact, if we calculate the curl of H40h1 taking into account (|41() . 
we get 

duj i J?-. 

— = (oj • grad) v + gradi x divT — Lcurl [divT ) . (54) 
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When the fluid is incompressible, the addendum grade x divT vanishes, and in 
a two-dimensional flow we have also (u> • grad) v = 0, but the vorticity is still 
altered by the term —icurl (iuT E ) , which is different from zero except when (i) 
there exists a scalar field B 3 x i — ► w (x) such that gradir = div T E , or (ii) there 
exists a second order tensor valued field B 3 x i — > A (x) e Horn (Vec, Vec) 

such that curl ( r T E ) = (curlA) T . So that, in time-dependent states, even in 
a two-dimensional flow of an incompressible perfect complex fluid the vorticity 
can be altered with consequent generation of topological transitions. 

Experiments confirming the theoretical predictions of the present paper in 
some special case could be useful toward the understanding of complex fluids. 
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